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Abstract 

Log-periodic oscillations have been found to decorate the usual power law behavior found to describe the 
approach to a critical point, when the continuous scale-invariance symmetry is partially broken into a discrete-scale 
invariance (DSI) symmetry. For Ising or Potts spins with ferromagnetic interactions on hierarchical systems, the 
relative magnitude of the log-periodic corrections are usually very small, of order 10^^. In growth processes (DLA), 
rupture, earthquake and financial crashes, they are found of the order of 10%. We offer a "technical" explanation 
for this 4-order-of-magnitude difference based on the property of the "regular function" g{x) embodying the effect 
of the microscopic degrees of freedom summed over in a renormalization group (RG) approach F{x) — g{x) + 
^~^F{'-^x) of an observable F as a function of a control parameter x. For systems for which the RG equation 
has not been derived, the previous equation can be understood as a Jackson g-integral, which is the natural tool 
for describing discrete scale invariance. We classify the 'Weierstrass-type" solutions of the RG into two classes 
characterized by the amplitudes An of the power law series expansion. These two classes are separated by a novel 
"critical" point. Growth processes (DLA), rupture, earthquake and financial crashes thus seem to be characterized 
by oscillatory or bounded regular microscopic functions that lead to a slow power law decay of An, giving strong 
log-periodic amplitudes. If in addition, the phases of An are ergodic and mixing, the observable presents self-affine 
non-differentiable properties. In contrast, the regular function of statistical physics models with "ferromagnetic"- 
type interactions at equibrium involves unbound logarithms of polynomials of the control variable that lead to a 
fast exponential decay of An giving weak log-periodic amplitudes and smoothed observables. These two classes 
of behavior might be traced back to the existence or abscence of "antiferromagnetic" or "dipolar"-type interactions 
which, when present, make the Green functions non-monotonous oscillatory and favor spatial modulated patterns. 
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1 Introduction 



The existence of log-periodic oscillatory corrections to the power laws associated with critical phenomena and, more 
generally, to observables of systems endowed with the scale-invariance symmetry has been recognized since the 1960s 
(see [jl|] for a recent review and references therein). The log-periodic oscillations result from a partial breakdown of the 
continuous scale-invariance symmetry into a discrete scale-invariance symmetry, as occurs for instance in hierarchical 
lattices. 

However, for one of the most studied class of models exhibiting these oscillations, i.e.. Potts model with ferro- 
magnetic interactions on hierarchical lattices, the relative magnitude of the log-periodic corrections are usually very 
small, of order 10^^ [^. In contrast, in growth processes (DLA) |^, 0], rupture earthquakes and financial 
crashes [^, ^J, they are found of amplitude of the order of 10%. 

Here, we propose an explanation for this puzzling observation of an 4-order-of-magnitude difference based on the 
nature of the microscopic interactions of the systems. Within a renormalization group (RG) approach, an observable 
at one scale can be related by a functional relation to the same observable at another scale, with the addition of 
the contribution of the degrees of freedom left-over by the procedure of decimation or of change of scale. This 
contribution is called the "regular part" of the renormalization group equation of the observable. For systems for 
which the RG equation has not been derived, the RG equation can be understood without reference to the RG as 
a Jackson g-integral [Q], which is the natural tool [1^, 11] for describing discrete scale invariance. Here, we do not 



discuss the mechanisms by which the continuous scale invariance symmetry is broken to give discrete scale invariance 
but rather present a phenomenological approach based on the functional RG/Jackson g-integral equation. 

Using the Mellin transform applied to the formal series solution of the renormalization group, we identify two 
broad classes of systems based on the nature of the decay with order n of the amplitudes An of the power law series 
expansion of the observable: 

1. systems with quasi -periodic "regular part" and/or with compact support have coefficients An decaying as a 
power law An ~ n^^, leading to strong log-periodic oscillatory amplitudes; if in addition, the phases of An are 
ergodic and mixing, the observable presents singular properties everywhere, similar to those of 'Weierstrass- 
type" functions. 

2. systems with non-periodic "regular part" with unbound support have An decaying as an exponential An ~ 
g-Kn q£ jj^gjj. oj-fjej- leading to exceedingly small log-periodic oscillatory amplitudes and regular smooth 
observables. 

We find families of "regular parts" which belong to both classes, with a "critical" transition from the first to the other 
as a parameter is varied. 

A known example of a system of the first class is the g-state Potts model with antiferromagnetic interactions 



[12, 13]. Another example is the statistics of closed-loop self-avoiding walks per site on a family of regular fractals 



with a discrete scale-invariant geometry such as the Sieirpinsky gasket [|14|]. A known example of the second class 
is the g-state Potts model with ferromagnetic interactions [^. Based on this knowledge and on the analysis of this 
paper which underlines the importance of the oscillatory property of the "regular part" or of its compact support to 
generate strong log-periodic structure, we conjecture that systems that possess "antiferromagnetic"-type or "dipolar"- 
type interactions that introduce modulated spatial patterns in the system belongs to the first class. This provides 
a possible explanation for the strong log-periodic oscillatory amplitudes observed empirically in growth processes, 
rupture, earthquake and financial crashes as we argue in the concluding section. 

The next section 2 introduces the renormalization group with a single control parameter, its formal solution with 
the presence of log-periodic corrections in the presence of discrete scale invariance. Section 3 uses the Mellin trans- 
form to resum the formal series solution of the renormalization group into a power law series. Section 4 presents the 
general classification within the two classes alluded to above in terms of the leading exponential or power law decay 
of the coefficients of this power law expansion. It examines the conditions under which the observable can develop 
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non-differentiable fractal properties similar to Weierstrass-type functions. A family of "regular parts" is introduced 
which exhibits a critical transition between the two classes. Section 5 presents many more examples of both classes. 
Section 6 concludes with a qualitative argument explaining why rupture and financial systems should belong to the 
first class. The table 1 offers a synthesis of the classification in terms of the decay of the coefficients An of the power 
series expansion of the observable for various choices of the "regular part" of the renormalization group. 



2 "Weierstrass-type functions" from discrete renormalization group equations 

Speaking about a material shape, a mathematical object or function, the symmetry of scale invariance refers to their 



invariance with respect to changes of scales of observation (see [|15|,[lg] for general introductions). In a nutshell, scale 
invariance simply means reproducing itself on different time or space scales. Specifically, an observable / which 
depends on a "control" parameter x is scale invariant under the arbitrary change x ^ 7a; if there is a number /i(7) 
such that 

fix) = -fi^x) . (1) 
A* 

Such scale invariance occurs for instance at the critical points t = tc of systems exhibiting a continuous phase 
transition. The renormalization group theory has been developed to provide an understanding of the emergence of the 
self-similar property ([I]) from a systematic scale change and spin decimation procedure [17]. 



Calling K the coupling (e.g. K = c^/^ for a spin model where J is interaction coefficient and T is the tempera- 
ture) and R the renormalization group map between two successive magnification steps, the free energy / per lattice 
site, bond, atom or element obeys the self-consistent equation: 

f{K)=giK) + -f[RiK)], (2) 

where g is a. regular part which is made of the free energy of the degrees of freedom summed over between two 
successive renormalizations, ^ > 1 is the ratio of the number of degrees of freedom between two successive renor- 
malizations. In general, this relationship (^ is an approximation whose validity requires the study of the impact of 
many-body interactions. When these higher-order interactions can be considered secondary as the scale of description 
increases (corresponding to so-called "irrelevant" operators), expression (^ becomes asymptotically exact at large 
scales. For perfectly self-similar problems, for instance for physical systems with nearest-neighbor interactions de- 
fined on regular geometrical fractals such as the Cantor set, the Sierpinsky Gasket, etc., or on regular hierarchical 
lattices, expression (^ is exact at all scales. 
It is solved recursively by 

00 

fiK) = Y.-9[R^''HK)], (3) 

n=0 ^ 

where i?^") is the n*^ iterate of the renormalization transformation. Around fixed points R{Kc) = Kc, the renormal- 
ization group map can be expanded up to first order in K — Kc as R{K) = ^{K — Kc). Posing x = K — Kc, we 
have = 7"x and the solution (^ becomes 

00 1 

In principle, (^ is only applicable sufficiently "close" to the critical point x = 0, that the higher-order terms in the 
expansion R{K) = j{K — Kc) can be neglected. The effect of nonlinear corrections terms for R{K) have been 
considered in [ [l2| , ^. 

The form (js]) or (Q) has not been derived rigorously for growth, rupture and other out-of-equilibrium processes 
alluded to above, even if there are various attempts to develop approximate RG descriptions on specific models of 
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these processes. It may thus seem a httle premature to use this discrete renormalization group description for these 
systems. Actually, expression (Q) can be obtained without any reference to a renormalization group approach: as 



soon as the system exhibits a discrete scale invariance, the natural tool is provided by ^-derivatives [11] from which 
it is seen that expression ^ is nothing but a Jackson q-integral of the function g{x), which constitutes the natural 
generalization of regular integrals for discretely self-similar systems [jlT]]. The way the Jackson g-integral is related 
to the free energy of a spin system on a hierarchical lattive was explained in [[T^]. 

In the mathematical literature, the function is called a Weierstrass-type function, to refer to the introduction by 

oo 

/v^ = V6"cosK7rx] , (5) 



K. Weierstrass of the function 1 19] 



n=0 



corresponding to the special case fi = l/b, j = a and g{x) = cos[7ra;]. To the surprise of mathematicians of the 19th 
century, Weierstrass showed that the function (^) is continuous but differentiable nowhere, provided 0<6<l,a>l 
and ab > 1 + ^tt. Note that, in the context of the renormalization group of critical phenomena, the condition 
a = 7 > 1 implies that the fixed point Kc is unstable. Hardy was able to improve later on the last bound and obtain 



that the Weierstrass function (g) is non-differentiable everywhere as soon as a6 > 1 [|20|]. In addition. Hardy showed 
that it satisfies the following Lipschitz condition (corresponding to self-affine scaling) for at > I, which is much 
more than just the statement of non-differentiability: 

fw{x + h) — fw{x) ~ \h\"^ , for all x where m = ln[l/6]/lna . (6) 

Note that for ab > 1, m < 1, expression shows that fw{x + h) — fw{x) ^ \h\ for h ^ 0. As a consequence, 
the ratio [fw{x + h) — fw{x)]/h has no limit which recovers the property of non-differentiability. Continuity is 
obvious from the fact that fw{x + h) — fw{x) — > as /i — > since m > 0. For the border case a = b discovered by 
Cellerier before 1 850, fw is not non-differentiable in a strict sense since it possesses infinite differential coefficients 
at an everywhere dense set of points [pl]. Richardson is credited with the first mention of the potential usefulness 



for the description of nature of the continuous everwhere non-differentiable Weierstrass function []22p. Shlesinger and 
co-workers [ p3| ] have previously noticed and studied the correspondence between and Weierstrass function. 

If one is interested in the non-regular (or non-analytic) behavior only close to the critical point x = 0, the regular 
part can be dropped and the analysis of ([T]) is sufficient. It is then easy to show that the most general solution of (|I|) is 
(see [|l]] and references therein) 

where 

m = i^, (8) 

In 7 

and P{y) is an arbitrary periodic function of its argument y of period 1. Its specification is actually determined by the 
regular part g{x) of the renormalization group equation, as shown for instance in the explicit solution (^. The scaling 
law f{x) ~ X™ implied by (0) is a special case of ^ obtained by putting x = and replacing /i by x in (^). 

The Laplace transform of f{x) defined by (^ also obeys a renormalization equation of the type (^. De- 

noting giiP) the Laplace transform of the regular part g{x), we have 

hm = i:^^9L[pm, (9) 

and 

fL{f3)=gL{P) + —fL(-) . (10) 
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The general solution of ( [T0| ) takes the same form as (^: 

where Phiv) is an arbitrary periodic function of its argument y of period 1. 

3 Reconstruction of "Weierstrass-type functions" from power series expansions 

Following we use the Mellin transform to obtain a power law series representation of the Weierstrass-type 

function The Mellin transform is defined as 



oo 

s-1 



f{s) ^ / x'-'f{x)dx . (12) 
Jo 

The Mellin transform ( |T2| ) provides a reconstruction of the infinite sum of the Weierstrass-type function as a sum 
of power law contributions, j4„3;~*", with "universal" complex exponents s„ determined only by properties of the 
hierarchical construction and not by the function g{x), with amplitudes An controlled by the form of the regular 
part g{x). These "non-universal" amplitudes in turn control the shape of the function f{x), its differentiability or 
non-differentiability as well as its self-affine (fractal) properties, as we shall describe in the sequel. 
The Mellin transform of reads 

m = 9{s) , (13) 

where g{s) is the Mellin transform of g{x). The inverse Mellin transformation of f{s), 

1 rc+ioo ^ 

fix) = — / f{s)x~'ds , (14) 

ZTTt Jc— zoo 

allows us to reconstruct /(x) as a new expansion in singular as well as regular powers of x in order to unravel its self- 
similar properties. Indeed, the usefulness of the Mellin transform is that power law behaviors spring out immediately 
from the poles of /(s), using Cauchy's theorem. 

In inverting the Mellin transform, we have two types of poles. The poles of the Mellin transform g of the analytical 
function g{x) occur in general at integer values and contribute only to the regular part fr{x) of /(x), as expected since 



g{x) is a regular contribution. The poles of the first term in the r.h.s. of (13) stem from the infinite sum over 

successive embeddings of scales and occur at s = s„ where 

27r 

Sn = —m + i- n , (15) 

m7 

and m is given by Their amplitude An is obtained by applying Cauchy's theorem and is given by the residues 

1- exp(-27rnz) g{s) 

^™ — 's T = i 9{s) = -j • (16) 

s^s„ — 1 In 7 In 7 

The resulting expression for f{x) is 

f{x) = fs{x)+fr{x), (17) 

where the singular part fg (x) is given by 

oo 

fs{x) = ^ An X-'- , (18) 
n=0 
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and 

An = ^^. (19) 

This approach is similar to the one developed in [ ^ ] for "fractal strings" r] (for instance, the complementary of the 
triadic Cantor set is a special fractal string). Their fractal properties are fully characterized by the introduction of the 
"geometric zeta function" Cri{s), which can be shown to be nothing but the Mellin transform of the measure defined 



on the fractal string (see |25] page 73). In particular, the poles of Cn{s) give the complex fractal dimensions of the 
fractal strings, similarly to the role played here by the complex exponents Sn defined by ([l5l). 

The regular part fr{x) of f{x) defined in (lO) is generated by the poles of g{s) if any, located at s = —n. 



n = 0, 1.... The residues of these poles give the coefficients Bn of the expansion of the regular part as follows: 

oo 

fr{x) = Y,BnX'^. (20) 



n=0 



4 Classification of "Weierstrass-type functions" 
4.1 Classification 



The representation ( |18| ) offers a classification of "Weierstrass-type functions" as follows. We will work in the class 
of g{x) (not covering of course all possible types of behavior of A{n)) where the coefficients An can be expressed as 
the product of an exponential decay by a power prefactor and a phase 

An = j^^ e-'^" e'^- , for large n , (21) 

where p, k > and Tpn are determined by the form of g{x) and the values of fi and 7. This class is broad enough to 
include many physically interesting shapes of g{x) as will be illustrated at length below. 

4.1.1 Justification of tlie classification 



The parameterization ( |2T| ) can be seen to result from very general theorems on the Mellin transform [26, 27]. Let us 
assume that the function g{x) defined for x > is continuous and satisfies the following conditions 

\g{x\ < cix"^ , < X < 1; \9{x\ < C2 , 1 < x < 00 , (22) 

where a > p. Then, its Mellin transform is a regular (differentiable) function inside the strip —a < Re{s) < —fi, One 
should also bear in mind that Re{s) < because of the constraints imposed by the very formulation of the problem. 
All functions that we shall consider below as examples belong to the class of continuous functions satisfying slightly 
more restricted conditions [^] such as (^) with a > and /3 = 0. As a consequence, their Mellin transform is 
regular for —a < Re{s) < 0. For instance, g{x) = cos{x) — 1 corresponds to a = 1 and /? = 0. The same conditions 
apply to log{l + x) and exp{—x) — 1. For the stretched exponential function exp{—x^) — 1 with h > 0), we have 
a = h and again /5 = 0. 

We are interested, in g{sn) = g{—m + ina;), particularly as n goes to infinity. The general condition which is 



usually imposed on this quantity in order to ensure the existence of its inverse Mellin transform is [ 27 ] 



g{—m + inu) ^0, as n — > +00 . (23) 

Again, An must be designed in such a way that it satisfies this condition automatically. 

Let us consider some simple but vivid examples, intended to illustrate how a power-law and exponential decay of 
An as a function of n emerges from simple functions satisfying the conditions stated above. We also note that, when 
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g{t) possesses discontinuities of the first kind, it still yields the dependence (p|) of An as a function of n. Maybe the 
simplest function leading to An with a power-law decay is 

g{x) =0, < X < 1 and g{x) = -1 , 1 < x < oo , (24) 

which leads to g{s) = ^, which is regular within the strip oo < Re{s) < 0. The corresponding An decays in 
amplitude as for large n. Different strip geometries lead to the same power-law decay of An, for instance 

g{x) = x°-, < x < 1 and g(x) = 0, 1 < x < oo , (25) 

with Mellin transform g{s) = (s + a)^^ with —a < Re{s) < 0. Let us also consider 

g{x) = , < X < 1 and ^(x) = -x" , 1 < x < oo , (26) 

which leads to a similar Mellin transform g{s) = (s + a)"^ but a different strip geometry Re{s) < —a. The slightly 
more complicated example of a continuous function composed of power laws 

g{x) = {b - ay^x"- , 0<x<l and ^(x) = (6 - a)"^x\ l<x<6, (27) 

leads to g{s) = (s + a)~^{s + b)^^ with —a < Re{s) < —b and the amplitude of An decaying as n^^. The analysis 
of these examples and of their Mellin transforms at s = s„ demonstrate that particulars of the strip geometry in the 
variable s are not important when one is concerned with the large n asymptotic behavior of g{sn)- The asymptotic 
power decay of An as a function of n can be dominated by an exponential decay, as we shall see in more details 
below. For instance, the continuous function formed by compounded power-laws g{x) = ^ x(l + x)^^ leads to 

g{s) = — csc(7rs), — 1 < Re{s) < 0, yielding An decaying as exp{—n) as n ^ +oo. 



Violations of the parameterization (pT|) regarding An occur when the conditions of the theorem |26] are changed, 
e.g. when the argument x is replaced by, say, ln(l/x) or when singularities are introduced to the function ^(x). This 
can be seen from equation (19) which shows that A^ is proportional to the Mellin transform of g{t) expressed at 
s = Sm = —m + inuj where 

UJ = ^. (28) 

m7 



1 f+oo 

An = — duG{u)e''^^\ (29) 



Posing u = In X, the Mellin transform becomes a Fourier transform 

1 

In 7 
where 

G('u) = e"(i-™)5(e") . (30) 
It is clear that, by a suitable choise of g{x), any dependence of An can be obtained. For instance, for 

G{u) = u-=^/2 g-a/« ^ (3^) 



we obtain An ~ e"^^"^" cos[\/2an], which exhibits an oscillatory stretched-exponential decay intermediate between 
the exponential (k > 0) and pure power law decay (k = 0) of (pT]). However, the choice ( |3T| ) corresponds to a rather 
special choice for 

(lnx)-3/2 e-'^/i'i^ 

In this case, g{x) +00 for x ^ 0, and this case is outside the domain of validity ( p2| ) of the theorem |26, 
Consider also the following example 



9ix) = TTz;^ . (32) 



gix) = .-y^e.J-'^] , (33) 
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leading to g{s) = exp(s^), valid for arbitrary s, which leads to An with amplitude decaying as exp{—n'^) as n — > oo. 
This example is also characterized by a pathological behavior for x — > +00 of g{x) which diverges faster than any 
power law. Another pathological example is 

g{x) = 1/2 7r"^/2 cos ((1/4) ln{l/xf - 7r/4 ) , (34) 

leading to g{s) = cos s^, valid for arbitrary s, which yields An with an amplitude growing as exp(n) as n ^ +00. 
This violates the condition on the Mellin transform given in p7|]. 



The existence of discontinuities of g{x), as one might expect from the theorem [ ]26| , |27| ], also violates the pa- 
rameterization (^ of An- Consider g{x) = — ^ j;^/^(l — x)~^ with an integrable singularity, which gives g{s) = 
tan(7rs), —1/2 < Re{s) < and An with an amplitude bounded from below by a constant as n — > +00. 



This absence of decay allows us to reject this type of function, since a decay of An is required by [g7|j, In con- 
strast, a logarithmic singularity, as for instance in g{x) = ^ log , is allowed. In this case, this gives g{s) = 

tanh(7rs), —1 < Re{s) < and the amplitude of g{sn) exhibits periodic modulations as n — > +00 as 
(sinh^(^7rm) + cos ^(^vrwn)^ {sm{Trojn) + i sinh(7rm,)). Another example with the logarithmic function 



(93) discussed below gives a power law decay with a logarithmic correction as shown by (p5\) due to the presence of 



the singularity. 

In conclusion, as long as the conditions of theorem [26] on Mellin transforms hold, the dependence of An as 
n —>■ +00 given by ( |2T| ) will hold as well. Violations of the theorem due to a change of variable or to the presence of 
simple poles lead either to a faster decay or to a non-decaying An- Allowing for logarithmic singularities within g{x) 
brings in logarithmic or oscillatory corrections to A^ as a function of n. 

4.1.2 Beyond the linear approxiation of the renormalization group map 

The asymptotic expansion (^Tj) uses the linear approximation i?^") (x) = j'^x that allows us to transform the general 
solution @ into the "Weierstrass-type function" @. As we said, (^ is only applicable sufficiently "close" to the 
critical point x = 0, such that the higher-order terms in the expansion R{x) = 7X can be neglected. The linear 
approximation of R^"'\x) = 7"x is bound however to become incorrect as n becomes very large, i.e., in the region 
determining the singular behavior. As discussed in [^, 24], the crucial property missed by the linear approximation is 



that /(x) is analytic only in a sector |arg x\ < 6 while we treated it as analytic in the cut plane |arg x| < vr. This im- 
plies that the exponential contribution e^'^'"- of the true asymptotic decay of the amplitudes of successive log-periodic 
harmonics is slower than found from the linear approximation, and goes as e~'^^"'. The angle 9 depends specifically 
on the flow map R{x) of the discrete renormalization group |^ and is generally of order 1. Our classification in two 
sets K = and k / is not modified by this subtlety. Here, we shall consider only the "Weierstrass-type functions" 
@ and will revisit the impact of nonlinear terms of the renormalization group map in a future communication. 

4.1.3 K>0: C°°-differentiabOity 

The general solution (0) remains true for any choice of the regular part ^(x) with the exponent m given by (||). 
This implies that there will always be an order of differentiation sufficiently large such that its becomes infinite at 



X = [12]. This is the crux of the argument on the existence of the singularity at x = 0. Here, we investigate the 



differentiability of /(x) for non-zero values of x, i.e., away from the unstable critical point x ^ 0. Expression (18) 
with (|l5j) provides a direct way for understanding the origin of the singular behavior at x ^ 0, as x"^ is in factor of 
an infinite sum of oscillatory terms with log-periodic oscillations condensing geometrically as x — > 0. 

When K > 0, the modulus of An decay exponentially fast to zero. Hence, /(x) is differentiable at all orders. This 
can be seen from the fact that 



00 



dx^ 



{-Sn)i-Sn - 1)...(-Sn - i + I) An X"^""^ (35) 
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is absolutely convergent for any order I of differentiation. Taking into account that | = x'^~ is independent 

of n and can be factorized, the nth term in the sum has an amplitude bounded by a constant times exp[— gn] since 
{—Sn){—Sn — Sn — ^ + 1) is bouudcd from above by a constant times n^. The sum is thus controlled by the 

exponentially fast decaying coefficients An and converges to well-defined values for any I. As a consequence of the 
exponential decay of the coefficients An, the log-periodic oscillations are extremely small. 

Another obvious way to ensure differentiability even when k = (see next section) is to truncate the number n of 
powers in the sum ( [T8| ) to a finite value: 

N 

f^s''\x) = J2AnX-'" , (36) 

n=0 

An example with iV = 1, 2, 3 is shown in figure 1 for the Weierstrass function (a = 7r/2 and p = m + \). For = 1, 
the real part /i^^ {x) is given by 



/P(x) = ao 
where oj is given by (El) and 



An=Q jn I H^n=l m /' i / \ , \ 

IH X H X cos[ll! m[x) + if) 

ao ao 



ao = — ^ , (37) 
/i- 1 



' Ini (An=i) 

(f = arctan ( ) + ktt, k = 0, ±1... (38) 

-flC [An=l ) 



This expression (37) is based on the singular part ( [l8| ) of the Mellin decomposition of the DSI equation It applies 
not only to the Weierstrass function but to any function of the form (^. Keeping only the first two terms recovers 
exactly the log-periodic formula introduced in the study of precursors of material failure [ p8| , |5|], of earthquakes 
precursors ||6|, |T]] and of precursors of financial crashes 

4.1.4 Critical behavior and non-differentiability 

Expression (p^) with ( |l5| ) and ( ^I| ) shows that fs{x) has the same differentiability properties as 

n=l 

Changing variable a; — > y = In x/ In 7, this reads 

+00 -. 

n=l 

With respect to the differentiability property, it is sufficient to study the real part of the infinite sum which reads 

^..,.,.,fa) = g °"''""J^'^"' . (41) 

71=1 

This expression allows us to recover some important results in the case where the phases ^„ are sufficiently 
random so that the numerators cos[27rny + ipn] take random uncorrelated signs with zero mean. Then, the sum 
^p,{^n}iy) truncated at n = T has the same convergence properties for T — > 00 as 

^m=f^. (42) 
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where dWt is the increment of the continuous white noise Brownian motion of zero mean and correlation function 
{dWtdWf) = 6{t — t)dt where 6 is the Dirac function. We get {X{T)) = and its variance is 

{[X{T)f)= / {dWtdW,)t-H'-r^ = -5-, (43) 



1 Ji 



1 t2p 



which is finite for T —>■ +oo if p > 1/2. This entails the convergence for p > 1/2 of the infinite series (^) for most 
phases tpn which are sufficiently ergodic and mixing. We thus expect that Kp^^^^j{y) and as a consequence fs{x) 
are continuous functions for p > 1/2. We can proceed similarly for studying their ^'s derivative. With respect to 



the convergence property, taking the ts derivative has the effect of changing p into p — £ in ([42[). We thus expect 
^p,{ipn} (y) ^^'^ ^ consequence fg (x) to be differentiable of order ^ for p > £ + 1/2. 
We conjecture the following. 

Conjecture on conditions for nondifferentiability from the singular power law expansion of "Weierstrass- 
type functions": provided that (1) k = and (2) the phases ipn are ergodic and, using a generalization of Hardy's 



condition ab > I for the Weierstrass function, the smallest order ^min of differentiation of fs{x) (defined by ( |T8| ) with 
exponents s„ given by ( p3| ) with (^) which does not exist is such that 

1 3 

- < P - Cin < - , (44) 

i.e., 

Cin = Int[p - ^] (45) 

is the integer part of p — |. In particular, with ergodic phases ipn of zero mean, the function /<j(x) is nondifferentiable 
forp < 3/2. 

It follows from Lebesgue's theorem on continuous functions of bounded variations that a non-differentiable func- 
tion is not a function of bounded variation. Therefore, a non-differentiable function is everywhere oscillating and 
the length of arc between any two points on the curve is infinite [^. This explains the observation below that 
the regular part g{x) must contain oscillations or must exhibit a compact support (so that it has a discrete Fourier 
series) in order for f(x) to be non-differentiable or for some of its derivatives to be non-differentiable. Actually, 
Weierstrass-type functions (Q) are believed to have the same Hausdorff dimension 2 — m as the Weierstrass function 
(^ for arbitrary regular part g{x), as long as it is a bounded almost periodic Lipschitz function of order (5 > m [|^]. 
The examples organized below in two classes illustrate and make precise this condition on g{x). We indeed find 
that non-differentiability occurs at a finite order of differentiation only for functions g{x) which are periodic or with 
compact support. 

It appears however that there is not yet a general understanding whether there exists a necessary and sufficient 
condition for the differentiability of a function on an interval. It is well-known that continuity is necessary for dif- 
ferentiability but is not sufficient as shown by the Weierstrass function and other examples above. The restriction 
of bounded variations has also proved insufficient: although a continuous function must possess a differential coef- 
ficient almost everywhere, yet there are examples of such functions which do not possess differential coefficients at 



unenumerable everywhere dense sets of points []21|]. 



4.2 General condition for k = 

Let us consider a regular function g{x) which is either periodic with period X or with compact support over the 
interval [0, X] and zero outside. It can then be expanded as a Fourier series 

g{x) = ^ + ^ [ofc cos(27rA:x/X) + bk sm{2TTkx/X)] , (46) 
^ k=i 
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where ao, ai, 61, a^, 6^, ... are arbitrary real numbers. 

The behavior of the coefficients An is controlled by the Mellin transform g{sn) of g{x) as shown by ([l9|). For 
g{x) periodic with zero mean, oq = and 



9{Sn) = 



+00 



+00 



fc=l 



fc=l 



(47) 



where cos(s) and sin(s) are the Mellin transform of cos x and sin x. Now, a general theorem on the Fourier series of 
periodic functions tells us that, if g{x) has continuous derivatives up to order r — 1 included and if the derivative of 
order k obeys the Dirichlet conditions, then the coefficients aj. and bj^ decay for large k as i.e., there is finite 

M' > M > such that M' > |afc| > M/k''+'^ and M' /k'^'^ > > M/k''^'^. If g{x) is discontinuous at a 
discrete set of points, this corresponds to taking r = in the previous formula. The Dirichlet conditions are: (i) g{x) 
is continuous or possess only a finite number of discontinuities; (ii) each point of discontinuity Xd is a discontinuity 
of the first kind, i.e., it is such that the limits to the left g{x x~^) and to the right g{x — > x~^) are finite; (iii) the 
interval [0, X] can be divided into a finite set of subintervals on each of which g{x) is monotonous. 
We can thus write 



+00 



+00 



k- 



smlSr 



' / t.r+1 — m 
k=l 



(48) 



The sum 



+00 

k=l 



{r+l—m+inuj) 



(49) 



is nothing else but the celebrated zeta-function ({y) of Riemann [ p5t |36| ], with the correspondence y = a + it, a = 
r + I - m, t = nuj. Ith known [||, ^ that \C{a + it)\ < {\t\ + l)V2-cr a < 0, where decreases like 
(27re)'^-i/2 for fj ^ — 00, and it does not satisfy a better estimate in this half -plane. For < cr < 1 (corresponding to 
< m < 1 andp = 0, C,{cj + it)\ < Kt^^^'^^^'^ ln(t)) uniformly for some constant K. However, we need the behavior 
of C(<7 + it) for o" = r + 1 — m > 0. It is obtained by using the relation ({s) = 2'*7r'*^^ sin ^r(l — s)C(l — s), 
which which can be separated into two parts that can be evaluated. Namely, sin(|s) r(l — s), re-casted in the 
variable z = 1 — s, takes the familiar form,cos(|z) T{z), which behaves for large n as n"''"'"™"^/^. The other term 
C{1 — s) = C{z) can be evaluated using the expression presented above for the zeta-function of an argument with 
negative real part (in our case for negative 1 — a and large n), ({z) < C {\n\ + i)^/^^^^^'^) = C{\n\ + i)i/2+»--m^ 
Therefore, the product of these two terms is of the order of C and the whole sum decay is slower than exponential. 

This shows that the sum (19) is of order 0{l/n^'^^^^'^^""^) and thus g{sn) is asymptotically a negative power of 
n for large n. This demonstrates that any periodic continuous function g{x) leads to a power law decay for An as a 
function of n. The continuity of g{x) implies that r > 1, ensuring that r + (1/2) — m > for m < 1. 

The same approach can be used for g(x) not periodic but defined on a compact support [0,X]. The discrete 
Fourier series expansion (46) still holds for x G [0, X] while g{x) = for x outside. A similar expression to ( p7| ) then 
holds in which oq 7^ in general and in which the Mellin transforms cos(s) and sin(s) are defined over the interval 
[0,X]. 



4.3 Bifurcation from wild to smooth "Weierstrass-type functions": an example using damped oscil- 
lators for the regular part of the renormalization group equation 

As a first example, let us consider the regular part g{x) of the renormalization group equation defined as 



9[x) 



■cos(a)x ^Qg^^ sin(a)) , with a £ 



0, 



(50) 
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The parameter a quantifies tlie relative strength of the oscillatory structure of g{x) versus its "damping": for a = t:/2, 
@ with (|^) recovers the initial function (|5|) introduced by Weierstrass with 6 = a = 7 and cos(7rx) replaced 
by cos(x); for a = 0, g{x) = exp[— x] has no oscillation anymore and corresponds to a pure exponential relaxation 



considered in |33] 



Plugging @ in (Q) gives 



oo 



where 



/(^) = E COS (7"^ sin(a)) , (51) 

n=0 ^ 

D = 2- m = 2- i^. (52) 
In 7 

The exponent D turns out to be equal to the fractal dimension of the the graph of the Weierstrass function obtained 
for a = 7r/2. Recall that the fractal dimension quantifies the self-similarity properties of scale invariant geometrical 
objects. Note that l<D<2asl<//<7 which is the condition of non-differentiability found by Hardy [|^] for 
the Weierstrass function. The graph of the Weierstrass function is thus more than a line but less than a plane. For 
a < 7r/2, /(x) is smooth and non-fractal {D = 1) and its graph has the complexity of the line. Actually, there are 
several fractal dimensions. It is known that the box counting (capacity, entropic, fractal, Minkowski) dimension and 



the packing dimensions of the Weierstrass function are all equal to D [37] given by (52) for a = 7r/2. It is conjectured 
but not proved that the Hausdorff fractal dimension of the graph of the Weierstrass function obtained for a = 7r/2 is 
also equal to D given by (^). It is known that the Hausdorff dimension of the graph of /(x) does not exceed D but 



there is no satisfactory condition to estimate its lower bound []34[]. 

Figure 2 shows the function ( pl| ) for a = tt/2 = 1.5708 (pure Weierstrass function: panel a), a = 0.9937r/2 = 
1.56 (panel b), a = 0.97r/2 = 1.414 (panel c) and and q = (panel d). 

The Mellin transform of g{x) defined by (^) for —1 < Re[s] = —m < (which is the interval of interest, as 
seen from (15)) is [38] 



g{s) = T{s) cos(a s) , (53) 

where T{s) is the Gamma function reducing to T{s) = (s — 1)! for integer arguments s. The additional term — 1/s 
disappears for <Re[s]. For values of the exponent m larger than 1, i.e., Re[s] < —1, additional correction terms 



should be added to ( p3[ ) [|38p. These additional terms only contribute to the power law dependence of the amplitudes 
An and not to the exponential. This problem is absent when the cosine in the definition of g{x) is replaced by the sine 
function. 

As we shall discuss below, the modification of g{x) into the modified function 

g;,j{x) = e-™^(") ^ COS {x sin(Q)) - 1 (54) 

gives g{s) = r(s) cos(a s) without the correction —1/s for —1 < Re[s] = —m < and leads to the so-called 
Mandelbrot-Weierstrass function. Similar "counter-term" should be introduced for stretched exponential and in simi- 
lar cases. They do not bring any extra contributions to the Mellin transform. 



The regular part fr{x) defined by (20) of f{x) defined in ( |T7| ) corresponding to g{x) defined by ( |50| ) is generated 
by the poles of T{s), located at s = — n, n = 0, 1..., since r(s) is analytic on the whole complex plane excluding 
these simple poles [^9|l. Using the expression Ress=-n r(s) = , we obtain its explicit form ( pO| ) with 

B{n)=^-^ ^cos(an). (55) 

n! fj, — ^"^ 

Note the particularly simple expression of the first term Bq = For |x| <^ 1, this constant term provides the only 
non-negligible contribution of the regular part fr{x) to f{x), whose behavior is completely controlled by the sum 
fs{x) of singular power laws. 
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The amplitudes An defined by (|T^ corresponding to g{x) defined by (5C) are 

r(Sn) COs(a Sn) 



In 7 



(56) 



The singular part fs{x), which is defined by (18) where the exponents s„ are given by (|T^), satisfies the exact scale- 
invariance equation (|l|). 

The asymptotic behavior of the amplitudes A„ given by (|5^ is 



Ania) 



n — > oo 



(57) 



with m = The angular log-frequency uj is defined by (|2^. In order to obtain (57), we have used the asymptotic 
dependence of the F— function asymptote for complex z 



T{z) ~ e 



{z-l/2) \nz-z 



Expression (^) is of the form ( pi] ) with p 
For a = 0, 

^n(O) ~ 



m + 



UJ 



\z\ > 1. 



a) and ipn = ^ n ln{uj n) 



n 



g— -^o; n UJ n ln(a; n) 



n 



oo 



(58) 



(59) 



As we have shown above, the fast exponential decay of A„(0) ensures the differentiability of f{x) at all orders. 
Actually, the fast decay of An{0) washes out any observable oscillatory structure from the function as seen on figure 
2d. However, there are very tiny log-periodic oscillations of amphtude less than 5 • 10^^ (see table 1) which are 
however unobservable at the scale of the plot of figure 2d. 

For a = 7r/2 (Weierstrass function), the exponential part disappears and 



An{7T/2) 



1 



n 



UJ n \n(u) n) 



n 



oo 



(60) 



This situation corresponds to the case p = m + ^ and k = and ipn = ^ n \n{uj n) in expression (^ of the 
classification of section 4J_ . The cancellation of the exponential term in A^. is due to the very peculiar compensation 



of the exponential decay of r(s„) by the exponential growth of cos(a s„) in (53), which occurs only for a = ^. 

The original Weierstrass function (Bh is thus seen as a very special "critical" or bifurcation point of the class of 
"Weierstrass-type functions" (Q) with (|50|). The analogy goes further as the expression (21) for the amplitudes An has 
the same structure as the correlation function of a system of spins where the order n in the sum (|18|) plays the role of 
the distance r between two spins. In this analogy, the "correlation length" is proportional tol/K~(| — q) 
diverges at the critical point a = f . 



,-1 



and 



4.4 Role of the phase: locaUzation and delocaUzation of singularities 



The phases tpn defined in (21) play an essential role in the construction of the self-affine nondifferentiable structure 



of the "Weierstrass-type functions". To stress this fact, let us consider several cases using different phases ijjn with 
the same absolute values \An\ of the amplitudes. This study parallels in a sense that of Berry and Lewis [41] and of 



Hunt [ 142[ ] but is distinct from it in a essential way as the phases considered here decorate the amplitudes An in ( |18[ ) 
of the power series expansion, rather than the phases in the cosine in (^. Actually, Berry and Lewis study a slight 
modification of the Weierstrass function (B|) defined as 



fwM = ^IP{1- cos[a"7ra;]) 



(61) 



n=0 
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proposed by Mandelbrot [43], which has the property of directly satisfying the "self-affine" property (|l|) with // = 1/6 
and 7 = 1. As discussed above, the choice (|^) for g{x), which gives (^) up to a sign, has the advantage of getting 
rid of the — 1/s correction in its Mellin transform ( |53| ) which makes thus more apparent and direct its self-similar 
properties. 



Hunt [ 42 1 is able to show that, by replacing the argument o^ttx of the cosine by a"7rx + 9n where 0„ are uncor- 



related random phases, the Hausdorff dimension of the phase-randomized Weierstrass function is D = 2 — m. 
4.4.1 Localization of singularities 

Let us first study the case where V'n is put equal to 0, i.e., we construct a phase-locked Weierstrass function as 

oo 

M^=J2\An{7T/2)\ X-'" , (62) 
n=0 

i.e., by constructing the singular part as the sum over power laws with amplitudes equal to the modulus of the ampli 

Os((71 

In 7 



tudes ( |56[ ) obtained for the Weierstrass function with a = tt/2, i.e., |A„(7r/2)| = | r(sn)coh((7r/2) s„) |^ without the 



phase. As a consequence, (60) is changed into 



A„(7r/2)| = C , forn ^ oo , (63) 



where C is a constant. 

Figure 3 shows the function fs{x) defined by (|62|) for m = 0.2 (panel a) and m = 0.65 (panel b). Rather than the 
familiar nondifferentiable self-affine corrugated structure of the Weierstrass function, fs{x) seems to be differentiable 
everywhere except for a discrete infinity of spikes at positions Xu, where u is an integer running from — oo to +oo, 
organized according to a geometric log-periodic structure. This discrete set of spikes decorates the leading singular 
behavior f{x) ~ for x — > of the general solution (^. Note that, in this case, the periodic function P of 
the general solution (^ is formed by the set of spikes geometrically converging to the origin. 

The spikes seem to diverge for m = 0.2 while they converge to a finite value for m = 0.65, as far as the numerical 



construction suggests. Appendix A examines some differentiability properties of (|62|). Appendix B shows that the 
functional shapes of the spikes for x — > = 1 /7" with u integer are given by 

Gs{x) . (64) 

|X-X„|2 

Thus, for < m < 1/2 (panel (a) of figure 3), the spikes correspond to a divergence of Gs{x) as x ^ x^. For 
1/2 < m < 3/2, Gs{x) goes to a finite value as x — > x^ but with an infinite slope (since < m — ^ < 1) according 
to Gs(x) constant — |x — x„|"^~2. 

Another example of "localization of singularities" is provided by the function 

oo 

/^(x) = J2 ^"""^ X-'" . (65) 



Figure 4 shows this function /s(x) defined by ( |65| ) with m = 0.2, to = 7.7. One can observe a log-periodic set of 
structures, each structure composed of log-periodic oscillations converging to singular points beyond which damped 
oscillation can be observed. Here, the phase Tpn = ^ ln{LO n) is not varying fast enough with n to scramble the 
complex power laws x~*", except at isolated points. 
Figure 5 shows the graph of 



oo 

, ^ ^ 1 



f^[x) = }_J n"'"-t e'"^ X-'" , (66) 

n=l 
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with m = 0.2, w = 7.7 and phases ipn = ^n. Again, the phase are not sufficiently random to make the function 
irregular, except at isolated points where the constructive interference of the phases lead to the isolated singularities. 



Note that both functions (p^) and ( pq ) can be analyzed with the method of Appendix B to obtain the functional 
form of the singularities. 

4.4.2 Mixing phases 

In contrast to the previous examples where the phases V'n are too regular, let us now take 

=ojn \a{uj n) (67) 



corresponding to the asymptotic dependence (pO|) of the amplitudes An of the Weierstrass function. The phases ( |67| ) 
are ergodic and mixing on the unit circle. The corresponding function is 

oo 



S[x) = """"^ " '■"('^ X-'- , (68) 

n=l 

which we call the "log-periodic Weierstrass" function to stress the fact that it is constructed by summing log-periodic 



power laws x (see for instance expression (37)) with amplitudes determined by the asymptotic behavior of the 



amplitudes of the power expansion of the Weierstrass function itself. The exponents s„ are again determined by (15) 
with m given by (|8|). By constructing (^), we are stripping off the Weierstrass function of its regular part and of 
all features that are unrelated to its fundamental nondifferentiability and self-affine properties. The definition of this 



"log-periodic Weierstrass" function ( |68| ) and its many generalizations studied below examplifies the novel construction 



developed here. S{x) exhibits the same non-differentiability as does the Weierstrass function. In all these cases, S{x) 
ondifferentiable since p = 
Consider the general case 



is nondifferentiable since p = m + ^ < |,in agreement with the conjecture of section ^ 



oo 



Si{x) = Y n"""^ exp (i V'W^ X-'" . (69) 
As other examples, let us now take 



n=l 

',(1) = ,,,^2 



i^'n> =u;n\ (70) 



and 



V'l^) = (J e"/"^ (71) 

and form the corresponding sums for i = 1 and 2. 'il^n\'ipn^ and ^pn^ cause similar irregular oscillations of cos[V'n^] 
between —1 and +1 as a function of n, allowing for a non-trivial and complex interactions of singularities whose 
amplitudes (most important as we have seen) exhibit a slow power-law decay. As the result Si{x) and 52 (x) exhibit 
very clear non-differentiable features shown in figures 9 and 10. 

Another very simple example of an ergodic phase is the quadratic rotator with irrational rotation number < i? < 

1: 

= + 27TRn . (72) 

The most irreg ular phase is obtained for R = g = = 0.61803398875... which is the golden mean whose 

main property is that it is the least-well approximated by a rational number. The corresponding "Golden-mean- 



log-periodic Weierstrass function" slf\x) defined by (|^) with ( |72[ ) is shown in figure 11. Other examples with 



R = 7r/4 = 0.785398163... and i? = 1/e = 0.367879441... lead to S^''^^\x) and S^^^''\x) shown in figures 12 and 
13. To each irrational number R corresponds an interesting "Weierstrass-type functions" whose delicately corrugated 
self-affine structure is encoded in the number-theoretical properties of its corresponding irrational number R. 

Note that if ( \J2\ ) is changed into ip^+i = V'r» + 27rii, the corresponding observable f{x) becomes smooth almost 



everywhere except at isolated points (also organized according to a geometrical series as described in appendix B). 
The ergodic but non-mixing properties of the linear rotation map does not scramble the phases sufficiently to create 
non-differentiability. 
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5 Illustration of the classification 



Different physical problems will be encoded by different regular parts g{x) quantifying the impact on the observable 
f{x) of the degrees of freedom summed over between two successive magnifications with the ratio A. For a given 
physical problem, we perform the Mellin transform of g{x). Then, together with the renormalization group structure 
@, its formal solution (||) and its expansion (^) close to the leading critical point x = leading to "Weierstrass-type 
functions", the classification (^Tj) of the previous section allows us to characterize the possible non-differentiability 
and scaling properties of the observable /(x). 



5.1 Other examples of the C°° differentiable family k> 

For a general statistical mechanics model, the regular part g{x) of the free energy has generally the form of the 
logarithm of a polynomial in x. Factorizing the polynomial, we do not lose generality by considering g{x) given by 

gix) = ln(l + x) , (73) 

for which 

g{s) = — ^ . (74) 

S Sm STT 

The poles of g occur for s = —n, n > 0, and contribute as already described only to the regular part of f{x). 
Therefore, the Mellin transform of /(x) is 

/(») = ^ -V^ ■ P« 

— Is sin(s7rj 

The regular part, determined by the poles at s = —n, n > 0, reads 

fr{x) = J2 B{m) X- , B{m) = ^ . (76) 

m=l ^ ' 

Note that constant term is absent. The singular part is 

oo 1 

/,(x) = Y.An x-^" , An = ^ . (77) 

In 7 sm(7r Sn) Sn 

The coefficients An converge to zero extremely fast for large n 

An ~ — e e , n — >oo, (/o) 

n 

where uj is given (28). This corresponds to p = 1, k = vr and ipn = —irm in the general classification (pT|). 
For the "lorenzian", 

g{x) = (l + x^)"' , (79) 
the coefficients An of the power law expansion ( f?^ of the singular part fs{x) are 

An = ^ ■ l \ ^ ' (80) 

2 m7 sm(7r/2 s„j 

with an asymptotic behavior given by 

An^ e-f'^'^ f n^oo. (81) 
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This corresponds to p = 0, k = ^io and ipn = general classification (pT|). The exponential decay rate 

^u> in this Lorentzian case is half that for the logarithm function (73). Both lead to C°° differentiable functions with 



extremely small amplitudes of log-periodic oscillations (see table 1). 
For so-called stretched exponential functions 



h>0 



(82) 



we obtain 



and 



An{h) 



r jSn/h) 
ln(7) h 



Anih) 



n h 



n 



oo 



(83) 



(84) 



Log-periodic stretched exponential function fs{x,h) and all its derivatives (on x) converge. This corresponds to 
p=^ + i, K = ^^ and V'n = ln{ujn) in the general classification (^Tj). Since k > 0, the corresponding 
singular function fs{x) is differentiable at all orders. However, a limit of non-differentiability at isolated points can 
be reached formally by taking the limit /i — > oo for which p ^ ^, k ^ 0^ and and ipn 0. Then, fs{x) exhibits 
the non-differentiability at points Xu verifying Xu = I/7" defined in (106) studied in section 4.4.1 . This is shown 
in figures 14 with the dependence of fs{x, h) on the parameter h. As h increases, g{x) becomes more and more 
localized close to the origin and fs{x) exhibits more and more pronounced steps. Formally, the limit /i — > 00 allows 
us to cross-over from the class k > to the class k = 0. 



5.2 Other examples of the "Weierstrass-type function" class k = 

We have noted above that a non-differentiable function is everywhere oscillating and the length of arc between any 
two points on the curve is infinite [^]. Its regular generator g{x) must thus contain oscillations or must exhibit at 
least compact support (so that it has a discrete Fourier series) in order for f{x) to be non-differentiable or for some of 
its derivatives to be non-differentiable. We illustrate this remark by several examples. 



5.2,1 Generalized periodic processes 

Let us consider the function generalizing the sinc-f unction by taking an arbitrary real exponent 5: 

.«) = 2^ . (85) 

The Weierstrass function is recovered for (5 = 0. The coefficients An of the expansion in power series of the singular 
part can be obtained by a simple shift of s in the expression obtained for the Weierstrass function: 

T{sn - S) sin ( iyr (s„ - 6)) 
An{5) = ^ ^ . (86) 

For m + 5 < 1, the log-periodic generalized sinc-function fs{x, 6) as well as its associated Weierstrass-type function 
are continuous but non-differentiable. 

The asymptotic behavior of the coefficient An is 

Anid) - n-'''~^~^/^e~' " , n^oo. (87) 



This corresponds top = m + 5 + l/2, k = and ipn = '^n ln(a;n) in the general classification (pTj). Figure 15 shows 
the generalized sinc-function fs{x, 6) for different values of 6 with for 5 = —0.1 (solid line), 6 = (dashed line) and 
6 = 0.1 (dotted line) for m = 0.4 and to = 7.7. 
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Another interesting case is the sine integral, g{x) = Si{x) = J^^ dv The coefficients of the power law 
series of the singular part are given by 

_ r(s„)sin(i7r s„) 



In 7 



(88) 



A m— 3/2 —i LO n ln(ajn)+^m 



n 



oo 



(89) 



with asymptotics 

corresponding to p = m + 3/2, k = and ^/;„ = u;nln(u;n) + in the general classification (^. All three 
functions /(x), and S{x) defined by (68) have a continuous but non-differentiable first derivative for m < 1. 

However, the delicate log-periodic conjugations are enhanced in the graph of S{x). 



5.2.2 Localized processes 

Let us now study functions g{x) with compact support such as 

g{x) = (1 - x^Y^'^ , < X < 1 , g{x) =0 for X > 1 with h>l and > 2 
The coefficients An of the power series expansion of the singular part Fs{x) are 

B (u, Sn/h) 



An{u,h) 



ln(7) h 



(90) 



(91) 



where B{x, y) = ^^^j^y-^ is the beta-function. Figure 17 shows the function /(x) obtained from the direct sum (Q). 
Figure 18 shows how the shape of the log-periodical structures steepen with increasing h, as the function g{x) evolves 
from a half-n shape to a the plateau g{Q < x < 1) = 1 and otherwise. The log-periodic geometrical series of 
plateaux and steps shown in figure 1 8 is reminiscent of the structures found for rupture [^, ||] and earthquakes 
^ precursors. 

The asymptotic behavior of ^4^(1^) is 



An{v, h) 



with z/ > 2, 



n 



oo . 



(92) 



corresponding to p = v, k = and V'n = ''i" in the general classification (^Tj). 

As the phases ^„ = tt are constant and their contribution can be factorized, the function /s(x, u) has a behavior 
similar to the function (|62|) analyzed in section 4.4.1. In particular, we recover the fact that the points Xu = I/7" 



given by (106) make the imaginary contribution of x vanish. As a consequence, they are the most singular points. 



An analysis similar to that presented in section 4.4. 1 can be performed. 
Another example corresponds to the logarithmic function 



g{x) = ln(l - x) , 



< X < 1 and ^(x) = for x > 1 , 



(93) 



with compact support. Figure 19 shows the corresponding Weierstrass function (m = 0.5, u; = 7.7, N = 47). The 
coefficients in the power expansion are given by 



Ar, 



1 



(94) 



where is the logarithmic derivative of the Gamma function. 
The asymptotic behavior of the coefficients An is 



Arn 



ln(n) 



n 



exp 



arctan 



2 In(a'n) 



vr 



, (n — > 00) 



(95) 
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corresponding to p = 1, k = and V'n = arctan ^ ^ in(t^") ^ jj^g general classification (pT]). The logarithmic 
"correction" to the power law 1/n comes from the singularity at x = 1. This example illustrates a possible cause for 
a deviation from the classification (^Tl). Such modification however does not change the qualitative picture as they 
correspond to the next sub-dominant correction to the power law contribution to An- 



The relative amplitudes of the two first power law terms are given by 
0.5,a; = 7.7). 

For the non-singular compact logarithmic regular part 



0.143, 



0.086, (m 



g{x) = ln(l + x) , < X < 1 and g{x) = for x > 1 , (96) 
figure 21 shows the corresponding Weierstrass-type function /(x). The coefficients in the power series expansion are 

1 



log(2)4^(l + f) + i^(i + f 



(97) 



where ^ is again the logarithmic derivative of the Gamma function. The asymptotic behavior of the coefficients An 
is 

1 r / / ti)n.\ \ 

(98) 



An '■ exp 

n 



I arctan I — ) ''^ ) 



[n 



oo 



corresponding to p = 1, k = and ipn = arctan (— ) + vr in the general classification (pTh. The relative amplitudes 



of the two first power law terms are given by 



0.016, 



^n = 



7.738 10-^ for m = 0.5, w = 7.7. Thus, 



even if the asymptotic decay is almost the same as for (93) up to the logarithmic correction, the logperiodic amplitudes 
of the leading terms are a factor of 10 smaller 



6 Discussion 

This paper has studied the solutions of the equation (^, which can be understood as a renormalization group equation 
with a single control parameter or more generally as the Jackson g-integral describing discrete scale invariant (DSI) 
systems. We have put the emphasis on the factors controlling the presence and amplitude of log-periodic corrections 
to the leading power law solution. We have used the Mellin transform to resum the formal series solution of the DSI 
equation into a power law series and have presented a general classification within two classes: 

1 . systems with quasi-periodic "regular part" and/or with compact support present strong log-periodic oscillatory 
amplitudes. 

2. systems with non-periodic "regular part" with unbound support have exceedingly small log-periodic oscillatory 
amplitudes and regular smooth observables. 

In systems for which the renormalization group equation has been explicited, systems of the first class are associated 
with "antiferromagnetic" interactions. Systems of the second class occur when the microscopic interactions are dom- 
inantly "ferromagnetic". These facts plus the observation that strong log-periodic oscillations have been reported in 
out-of-equilibrium growth processes, in rupture, in earthquakes and in finance lead us to propose that strong logperiod- 
icity is associated with the presence of interactions with competing signs and long-range dipolar and multipolar form, 
favoring different "signs" and modulated patterns. Dipolar and/or antiferromagnetic interactions are well-known to 



lead to a rich behavior with striped patterns with alternating up and down spins (see [ ]450 and references therein) which 
are robust in the presence of disorder such as in the dipolar random-field Ising model [^, 47]. We conjecture that these 
spatial organization translates into a quasi-periodic "regular part" (or with compact support) of the renormalization 
group description. 

The existence of "antiferromagnetic" or "dipolar" interactions leading to coexistence of different "signs" and of 
modulated structures is documented in rupture and finance. In rupture phenomena, cracking or local breaking leads 
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to a redistribution of stress according to the elastic "dipolar" Green function which is non-monotonous and exhibit 
oscillations as a function of distance with local enhancement as well as screening alternating as a function of angular 
position. Similarly, an earthquake (corresponding technically to a "double-couple" force source within the elastic 
crust) redistributes stress according to a clover-leaf pattern with alternating + and — signs of the stress as a function 
of angular position. This superficial analogy has been shown to be quite precise in a mapping between the random 
fuse model and a dipolar magnet in which the state of the network damage is associated with a metastable spin 



configuration [)48[]. 

Similarly, the complex price time series observed in financial stock markets can be seen to result from the competi- 
tion between "value-searching" investors [^, 5C ] who track the fundamental value of the stocks and "trend-following" 



traders (see [51] and references therein) who imitate the crowd thus developing herding behavior. Roughly speaking, 
the former (later) type of traders sell (buy) when the price increases above the fundamental value and sell (buy) when 
the reverse occurs. Thus, "value-searching" investors follow interactions similar to "antiferromagnetic" coupling 
while the herding behavior resulting from the action of "trend-following" investors can be captured quantitatively by 
"ferromagnetic"-type interactions. 

We believe that the existence of non-monotonous interactions with competing signs may be a fundamental mech- 
anism of strong log-periodicity. Our analysis of its impact on the "regular part" of the DSI equation suggests that the 
"antiferromagnetic" and/or "dipolar" interactions are crucial ingredients at the origin of fractal patterns in nature. 

However, determining the observables that take the place of the equilibrium free energy for growth models which 
could be obtained recursively (via a renormalisation group transformation) - or equivalently, the identification of the 
meaning of the function g{x), remains an unsolved problem. Correlatively, the mechanisms in growth models that 
select the log-period of the observed log-periodic oscillations in the absence of a hierarchical lattice have not been 
fully elucidated (see however [Q] for a specific needle growth model and ^ for a review). In a way, our paper 
has attempted to provide an explanation of why the log-periodic oscillations may have a strong amplitude, without 
explaining really their origin! As a consequence, the program ahead of us is i) to identify clearly the mechanism(s) 
underlying the emergence of DSI, ii) derive the equivalent renormalization group equation. Only then, shall we be 
able to substantiate convincingly our conjectured mechanism in terms of "anti-ferromagnetic" interactions. 

Acknowledgments: We are grateful to A. Erzan for a discussion on Jackson's integral and for supplying the 
correponding references and to the referees' suggestions that helped us improve the manuscript. 
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Appendix A: Differentiability properties of the "localization of singularities" 



Using the asymptotic expression (|63|), we can write 



oo -| / 1 \ 

Re[;^] = G{x) +x"' y COS ( 2^nf^ ) 



(99) 



where G{x) = E"Lo l^n(vr/2) 



is a regular function. Re[/s(x)] denotes the real part of fs{x) and we have 



used (p^. The second term of the r.h.s. of (99), which can be called the singular part of Re[fs{x)] and is denoted 
Gs{x), is a sum starting at an index rir which is taken sufficiently large such that the asymptotic expression ( |63| ) holds 
to within any desired degree of accuracy. 

The singular part Gs{x) has the same analytical behavior as the function 



°° 1 1 

V — cos(ny) , p = m+- 
~ nP 1 



(100) 



n=l 



where y = 111 hi x/ In 7. This function is a special case of |^,^| (y) defined by (^) for -(/^n = 0. 

This function Kp{y) has been studied in the literature for special cases. When the real part of p is larger than 
1, the sum is absolutely convergent for all x. Restricting our attention to real exponents p, the series Kp{y = vr) = 
— Y^^=i ^~7^~> which corresponds toy = {2i + 1)tt where i is an arbitrary integer, is convergent for all positive p 
to Kp{y = vr) = (1 — 2^~P)C(p) [pO|], where ({p) is the Riemann zeta function. Obviously, Kp{y = 2tt) is infinite 
for p < 1 and we show below that Kp{y — > 2it) has a power law singularity. For p > 1, Kp{y) can be expressed as 



4r(p) 



sec 



Tip 



(101) 



where C(s, v) = X]£o[^ + ^l ** generalized Riemann zeta function 

For p > and except for the special value y = modulus 27r, Kp{y) is finite and differentiable. To see this, let us 
consider rational values of y/27r = r/q with q > 2, where r/q is the irreducible representation of the rational ?//27r. 
We can rearrange the series in (|100D into q sub-series as follows: 



Kp{y) 



/ 1 



+ E 



,j27r 2{r/q) 



00 / ^ 



+ 2 



^ + E 



J2n {r/q) 



+ E' 



i27r j(r/g) 



P 



k=0 



{kq + j. 



{kq-2)P 

9/00 ^i2nj(r/q)' 



(102) 



We expand l/(% + j)^' = (^"^1 " i^J + + 



and get 



i^p(y/27r = r^) = Re£ 



fc=0 



^ gi27r i(r/g) _ P_^j ^i2-K j(r/q) ^ PjP +^) ^ -2 ^i27r j(r/g) 



2(M 



2 



+ ... 



j27r i(r/g) 



(103) 

(1 - w^g)/{l -w) = 0. 

in ( |103| ) is identically zero. The other sums are non zero and finite. We thus get 



CalUng Wq the qth root of 1, i.e., Wg = e*^'^/'^, we have then Wg + Wg + ... + ^ + 
Hence, the first sum X]?=i e 



Kp{y/27T = r/q) = £ 



k=0 



Ck{r,q) 



(104) 
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where 



Ck{r, q) 



Re 



j{r/q) 



p{p+l){p + 2) .3 



(105) 



is bounded from above as A; — > +00. The expression (104) shows that Kp{y/2'K = r/q) is finite for any p > 0. Now, 
since rational numbers are dense among real numbers, i.e., any irrational number can be approached arbitrarily close 
by a rational number, by the condition of continuity, Kp{y) is finite everywhere, except for y/27r = r/q with q = 1. 
Differentiating the expression ( |100| ) gives the series X^^i ~T=^ sin(ny). By the same reasoning leading to (104), 



this derivative is bounded from above by a constant times Y^^=i which is convergent for p > \. This shows that 
Kp{y) is differentiable for p > 1, and thus Re[fs{x)] is differentiable for m > 1/2. This approach is not powerful 
enough however to treat the case m < 1/2. 
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APPENDIX B: Functional shape of the denumerable set of discrete singularities re- 
sulting from the "localization of singularities" 



We now examine the special case y/27r = r/q with g = 1. From the expression (99), the values Xu, which are such 
that In Xu/ In 7 is an integer —u, i.e., 

Xu = 1/7" , (106) 
make all the cosine terms in the infinite sum in phase and equal to 1. Thus, 

00 

G.(x„)=x™ ^ (107) 

n=nr 1 1' ^ 

which diverges for m, < 1/2. At the border case m = 1/2, the divergence is logarithmic. Similarly, dGs/dx\x=xu 
diverges for m < 3/2 as an additional power of n is brought to each term in the sum by taking the derivative. This 
and expression ( |106| ) explain the graphs of Figure 3. 

The functional shapes of the spikes for x —>■ Xu can be determined as follows. For x x„, cos (21:11^^^ = 

cos (tI^) + 0{e'^), where e = {x — Xu)/xu and 0{e'^) represents a term proportional to e^. Let us now construct 
and compare Gs{e) and Gs{Xe), where A is an arbitrary number. Up to first order in e, we have 



00 



cos 



I 2'Kn\e 



G.(A6)..x- ' ■ (108) 

n=nr IT' ^ 

Posing n' = Int(nA), Gs{\e) can be rewritten 

a(A6)^x™ V ^^i^^ = x-A-i V ^ ' . (109) 

n'=lnt[nrX] " n'=lnt[nrX\ " ^ 



Note the presence of the additional multiplicative term ^^-^^ — in the sum of ( |109| ). The numerator A"* 2 stems from 

replacing n by n' = Int(nA) in 1/ n™+ 2 . The other factor 1/A is the "Jacobian" of the change from n to n' = Int(nA). 
Expression (|109[) can then be rewritten 



1 



G,(Ae) « A'"- 2 Gs{e) + Hr{e) , (110) 
where 

Int[n,.A]-l ( 27rn'e \ 

H{e) = X"'--2 y (111) 



n=nr 



n 



is a regular function of e. The singular part of Gs{e) is solution of Gs(Ae) A*" 2 Gs{e), i.e., 

Gs{e) ~ e"-^ . (112) 

This confirms that, for < m < 1/2 (panel (a) of figure 3), the spikes correspond to a divergence of Gs{x) 
according to Gs{x) ~ l/\x — x„| 2~"^. For 1/2 < m < 3/2, Gs{x) goes to a finite value as x ^ Xu but with an 
infinite slope (since < m — i < 1) according to Gs{x) ~ constant — \x — Xu|™~2 . The borderline case m = 1/2 
can actually be summed exactly as Kp=i{y) = — ^ In (2[1 — cos y]) [p^]. When y ^ modulo 2-7T, Kp=i{y) diverges 

as In i and thus Gs(e) diverges as Gs(e) ~ In 
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Table 1: Synthesis of the different classes of Weierstrass-type functions according to the general classification 

^ ^ e'^"" e^'^" of the expansion ( |T^ ) in terms of a series of power laws x~^". The parameters p, k > and ipn 
are determined by the form of g{x) and the values of fi and 7. All numerical values given in this table correspond to 
m = 0.5, Lo = 7.7 corresponding to 7 = 2.26 and fi = ^ = 1.5. The last two columns quantify the amplitude of 
the log-periodic oscillations with respect to the leading real power law. (*) c = cos a and s = sin a. 
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FIGURE CAPTIONS 



FIGURE 1: Power law expansion part fs given by (^) for the Weierstrass function (^), with = 1 (soHd), 
N = 2 (dash), = 3 (dot) oscillatory terms, respectively. Here, m = 0.25, u; = 6.3 corresponding to 7 = 2.7 and 
fi = 1.28. As the number of complex exponents increases, the number of the oscillations increase. 

FIGURE 2: Quasi-Weierstrass function for (a) a = f , (b) a = 0.9937r/2 = 1.56, (c) a = 0.97r/2 = 1.414 and 
(d) a = 0, for m = 0.25, to = 7.7, using = 32 terms to estimate the sums (^Ij). Increasing N does not change the 
results 



FIGURE 3: Panels (a) and (b) show fs{x) defined by ( p2| ) with zero phase ^„ = 0, for m = 0.2 and m = 0.65 
respectively, with the same lo = 7.7, constructed by truncating the sum at the = 29 term. The tiny regular 
oscillations result from the truncation to a finite A^ and slowly vanish when A^ ^ 00. They are thus spurious finite- 
size effects. 



FIGURE 4: Graph of /^(x) defined by (|5|) with V'n = w ln(a;n), with m = 0.2, uj = 7.7 and using A^ = 1000 
terms in the sum. 



FIGURE 5: Graph of shows fs{x) defined by (|6§) with Vn = ^n, with m = 0.2, oj = 7.7, with A^ = 1000 terms 
in the sum. 

',(1) 



FIGURE 6: Graph of Si{x) defined by (|6^ for the phases ipk ' defined by (^ with m = 0.25, lo = 7.7, N 
1000. 



FIGURE 7 Graph of S2{x) defined by (|69|) for the phases V'l^^ defined by ([7l| with m = 0.5, oj = 8, N = 200 



FIGURE 8: "Golden-mean-log-periodic Weierstrass function" S^^^ (x) defined by ( |69| ) with (^) for m = 0.5, w 
7.7, N = 500. 



FIGURE 9: "7r/4-log-periodic Weierstrass function" S^r^^^^ defined by (||) with (^) and i? = 7r/4 for m = 
0.5, = 7.7, N = 500. 

FIGURE 10: "e-log-periodic Weierstrass function" S^^^ defined by (||) with (^) and = e = 2.718... for 
m = 0.5, w = 7.7, A^ = 500. 

FIGURE 11: Singular part fs{x) of the Weierstrass-like function for the regular function g{x) equal to the 



stretched exponential (82) for /i = 5 (solid line), h = W (dashed line), h = 20 (dotted line), h = 50 (dashed- 



dotted line) and h = 100 (dashed-dot-dotted Une), for m = 0.4, lo = 7.7, N = 22. 



FIGURE 12: Singular part fs{x, 5) of the Weierstrass-like function for the regular function g{x) given by (85) for 
6 = -0.1 (solid line), 5 = (dashed Une) and 5 = 0.1 (dotted Une) for m = 0.4, uj = 7.7, N = 27). 



FIGURE 13: Weiersti"ass-type function f{x) for compact g{x) given by ( |90| ) with v = 2,h = 2fov m = 0.5, cj = 
7.7, A^ = 47). 

FIGURE 14: Evolution of the "singular part" fs{x) corresponding to the compact regular part g{x) (^) forv = 2 
with increasing abruptness of g{x) quantified by the exponent h: h = 2 (solid), h = 5 (dash), h = 10 (dot), h = 20 
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(dash dot), h = 50 (dash dot dot), h = 100 (short dash), for m = 0.5, io = 7.7, N = 47). 



FIGURE 15: Weierstrass-type function /(x) corresponding to the regular part g{x) defined by (?3) with compact 
support, with m = 0.5, lo = 7.7, N = 47. 



FIGURE 16: Weierstrass-type function /(x) corresponding to the regular part g{x) defined by ( p6| ) with compact 
support, with m = 0.5, lo = 7.7, = 47. 
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Fig. 13 

0.0 0.2 0.4 0.6 0.8 1.0 




0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 15 




0.0 0.2 0.4 0.6 0.8 1.0 

X 



